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Abstract 


A study of mode localization in mistuned bladed disks is performed using transfer matrices. The transfer 
matrix approach yields the free response of a general, mono-coupled, perfectly cyclic assembly in closed 
form. A mistuned structure is represented by random transfer matrices, and the expansion of these matrices 
in terms of die small mistuning parameter leads to the definition of a measure of sensitivity to mistiming. An 
approximation of the localization factor, the spatially averaged rate of exponential attenuation per blade-disk 
sector, is obtained through perturbation techniques in the limits of high and low sensitivity. The methodology 
is applied to a common model of a bladed disk and the results verified by Monte Carlo simulations. The 
easily calculated sensitivity measure may prove to be a valuable design tool due to its system-independent 
quantification of mistuning effects such as mode localization. 
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1. Introduction 

A perfectly periodic structure consists of a chain of identical elements connected to one another in an 
identical manner. If the last element of the chain is connected to the first one, the structure is said to have 
cyclic symmetry. When analyzing the dynamics of bladed-disk assemblies one frequently postulates perfect 
cyclic symmetry. This is equivalent to assuming that the bladed disk is tuned, that is, that all blades are 
strictly identical and are identically mounted and uniformly spaced on a homogeneous, symmetrical disk. 

All cyclic structures share the same set of modes of vibration, called constant interblade phase angle 
modes 1 . The vibration modes are so named, because in a normal mode motion all blades vibrate with the 
same amplitude and the difference in phase between the motion of adjacent blades is constant throughout the 
blade assembly. This assumption of perfect cyclic symmetry simplifies drastically the vibration analysis of 
bladed disks. Instead of analyzing the structure as a whole, the equations of motion may be uncoupled and 
the size of the problem reduced to that of a single blade or that of a blade-disk sector. However, designers 
are aware of inherent differences among rotor blades, due to material and manufacturing tolerances as well 
as in-service degradation, a phenomenon referred to as rotor mistiming. While it is widely accepted that 
mistiming causes a reduction in the risk of flutter instability [1,2], it has also been established that mistuning 
has a negative effect that could outweigh this benefit, by increasing forced response amplitudes: Thus 
mi stiming may alter the results of a tuned analysis drastically[3,4,6] 2 . 

One of the most important problems that plague turbomachinery rotors is the existence of rogue blades 
— lone blades that exhibit unexpected fatigue failure. Recently it has been suggested that rotor mistuning 
might be the cause of rogue blades, through a phenomenon called normal mode localization, whereby 
vibrations are confined to a few blades of the assembly [2-4,6-9]. 

The phenomenon of vibration localization may be expected to occur in any nearly periodic structure for 
which perfect periodicity is prevented by small irregularities. Localization, like damping, manifests itself 
as a spatial decay of the vibration amplitude along the structure, but for vastly different reasons. In the case 
of damping, energy is dissipated as vibrations are transmitted through the system, whereas in the case of 
loc alizati on, the energy is merely confined to a small geometric region within the structure. Localization 
occurs because waves propagating away from the energy source are reflected by the boundary between the 
slightly different subsystems making up the nearly periodic structure. The resulting confinement of energy 
may lead to much higher amplitudes locally than would be predicted if perfect periodicity were assumed, 
with possibly disastrous effects, for example in turbomachinery. The localization phenomenon has recently 
received wide attention in the literature and it has been shown to occur in various types of nearly periodic 
structures, namely blade assemblies [2-4, 6-9], multi-span structures [10-14], and some large space structures 
[15]. 

The manner in which the substructures, or bays, which make up a (nearly) periodic structure are 
interconnected plays a major role in its dynamics and in the occurrence of mode localization. It is only 
through these connections that vibrational energy is passed between substructures. If neighboring bays are 
connected through several degrees of freedom the periodic structure is said to be a multi-coupled or a multi- 
wave structure. Conversely, if neighboring bays are only connected through a single degree of freedom, the 
system is mono-coupled and it carries a single pair of left- and right-traveling waves [16]. 

Each bay in a (nearly) periodic system is, in general, a multi degree-of-freedom substructure which may 
or ma y not be connected to ground. In the work of Mead [16], a bay with at least one immovable point within 
is called “positive-definite” (for example, a multi-span beam on simple supports), whereas a “semi-definite” 
bay has no such immovable points (for example, a taut string with masses at regular intervals). In the case of 
a turbine rotor, each bay in the chain is a model of one blade and of an appropriate sector of the disk to which 
the blade is attached. These blade-disk elements are “positive-definite ” since they have fixed points in the 
form of either an unmodeled (rigid) part of the disk to which the elements are attached or their common 


1 In the periodic structure literature these inodes are often referred to as extended, or global modes. We shall make use of the latter terminology here. 

2 The numerous studies of the effect of mistuning on the dynamics of blade assemblies are reviewed in a survey paper by Srinivasan [5]. 
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axis of rotation. Coupling between blades is due to structural coupling through the disk and aerodynamic 
coupling through the fluid. Blades may also be connected through a shroud. Hence a blade assembly is, in 
general, a complex multi-coupled (nearly) cyclic structure. 

The dynamics of a tuned, multi-coupled blade assembly are modeled by matrices that have a circulant 
structure and are generally fully populated. Recall that a matrix is circulant if every line is a permutation of 
the previous line by a shift to the right by one element [17]. If several degrees of freedom are considered 
for each blade, the system matrices become block-circulant. The coupling between blades is generally 
strongest between neighboring blades, and becomes weaker as the blades are further apart. This leads 
to block-circulant matrices that are dominated by blocks close to the diagonal (and corresponding blocks 
near the top right and bottom left comers). The well-developed theory of circulant matrices [17] provides 
a valuable tool for the analysis of tuned cyclic systems. However, when cyclicity-breaking mistuning is 
introduced, the circulant structure disappears and another means of solution is desired. 

If the strength of coupling between two blades in an assembly decreases rapidly with decreasing 
proximity, a blade might be considered to influence only a few of its close neighbors, by ignoring the effect it 
has on the rest of the blades. To capitalize on this assumption, a transfer matrix approach can be introduced. 
For example, an assembly which features only coupling between adjacent blades can be modeled efficiently 
with transfer matrices of dimension twice the number of coupling coordinates [16]. Note that the number of 
degrees of freedom through which a blade is coupled to its neighbors may be much smaller than the actual 
number of degrees of freedom accounted for in the blade itself. With the transfer matrix approach the cyclic 
nature of the system is taken into account by realizing that the state vector describing the behavior of each 
blade in an AT-blade system is, at any given time, periodic with period N/n, where n is some integer. 

The present work centers on the study of the localization phenomenon in a general, mono-coupled, 
mistimed cyclic assembly. The primary goal of the study is the development of a measure of sensitivity 
to mistuning. Such a measure would ideally identify regions of high sensitivity in the space of the design 
parameters, thereby predicting the risk of the occurrence of strong localization in various classes of bladed- 
disk assemblies. The work is based entirely on a transfer matrix formulation of a mono-coupled assembly. 

The paper is organized as follows. In Section 2 the free dynamics of a general, mono-coupled, nearly 
cyclic structure are formulated with transfer matrices. Blade mistuning is introduced in the form of small, 
random parameters. For demonstrative purposes a blade assembly with two degrees of freedom per bay is 
introduced and used as an illustrative example for the remainder of the paper. Sections 3 and 4 contain, 
respectively, a study of the propagation of harmonic waves in, and of the natural modes of the tuned structure. 
In Section 5 we motivate the need for a sensitivity measure by illustrating the drastic effect that mistuning can 
have on the normal modes. Section 6 opens with a discussion of the effects of mistuning on the propagation 
of waves. The localization factor is subsequently defined as a means of quantifying localization. A measure 
of sensitivity is then suggested and analytical approximations of the localization factor are developed in terms 
of this sensitivity measure. The results are verified for the example system by Monte Carlo simulations. The 
main conclusions reached are given in Section 7. 

The primary contribution of the paper is a method for estimating sensitivity to mistuning. This is 
accomplished through the definition of a sensitivity measure, whose inexpensive calculation is based on 
the model of a single bay in a mono-coupled assembly. Furthermore, a quantitative, universal relationship 
between the sensitivity measure and the localization factor verifies the appropriateness of the sensitivity 
measure. 
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Figure 1 A general N - bay nearly cyclic assembly with P degrees of freedom per bay, one of which 
couples adjacent bays. 


2. Transfer Matrices for Nearly Periodic, Mono-Coupled Structures. 

The transfer matrix modeling of (nearly) periodic structures undergoing harmonic motion requires the 
definition of a state vector. For the discussion we introduce the terms bay and interface. A bay is exactly 
one spatial period in the periodic structure, such that only bays that are nearest neighbors are coupled. An 
interface is a point on the chain which separates two bays. To a state vector for a bay corresponds a transfer 
matrix that relates the states at two consecutive interfaces, or two consecutive bays. The dimension of the 
state vector must be twice the number of coupling coordinates at each interface [16]. A state vector is most 
commonly defined as the displacements of the coupling coordinates at an interface and the associated forces. 
For a mono-coupled structure (see Fig. 1) a two-by-two transfer matrix would relate the deflection and force 

at adjacent interfaces, thus: 

[f] ( = To 

Where the subscript 0 implies a tuned system. Alternatively we could define a state vector for a bay as the 
deflections of the coupling coordinates at both ends of the bay (see Fig. 1 ). Then 


i » 1, . . . , N. 


( 1 ) 


J i — 1 


Qi+l 
?. . 



t « 1, ... ,1V 


( 2 ) 


relates the states of two adjacent bays. This paper is based on the latter approach. The representation 
in Eq. (2) does not account for motion-independent external forces acting on the system and hence only 
the free dynamics are considered. The formulation of Eq. (2) requires two equations relating the coupling 
coordinates q i+ u q, and 3 . -1 ■ One is the equation of motion taken at interface i : 

-q i+ 1 + /3o(w)g» - Qi-I = 0 , fiM € JR. (3) 


Note that Eq. (3) is symmetric in g,+i and 1 , which requires the structure to be symmetric with respect 
to clockwise and counterclockwise numbering. This is the case only in the absence of aerodynamic and 
Coriolis forces. More on that later. The symmetry discussed here is different from Mead’s [16] definition 
of symmetry of individual bays, and Eq. (3) holds for either symmetric or unsymmetric bays. Equation (3) 
does not account for dissipation since /3 0 (u;) is assumed to be real valued. The second equation in Eq. (2) 
is simply the identity g, = g,. Hence, each bay of a perfectly cyclic (tuned) mono-coupled structure, is 
described by the same transfer matrix representation 



p 0 (v) e JR. 


(4) 


where Eq. (4) accounts for neither damping nor aerodynamic effects. In general each substructure possesses 
P degrees of freedom (see Fig. 1) and the generalized coordinates of the »th substructure are related to the 
coupling coordinate g t through P - 1 equations of motion local to the bay. However, note that the coupling 
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degree of freedom has no special significance, and in fact we may choose any of the P degrees of freedom 
as our reference coordinate q x . The choice must of course be the same for all bays. 

When mistuning is introduced into the above cyclic system its periodicity is broken in one of two ways. 
First, the mistuning may be caused by a parameter which only appears in relation to the z th interface, e.g., the 
mistuning of parameters of individual blades. When this is the case the symmetry of Eq. (3) is unaffected 
and Eq. (4) is replaced by 



where is the small deviation (order e or smaller) of the parameter from its average value, defining the 
mistuning for the I th bay. This is a random variable of mean zero. In the notation, the frequency dependence 
of (3 has been dropped, for clarity. 

Alternatively, the mistuning may be caused by a parameter which appears at both interfaces, z and 
(z + 1), e.g., the stiffness of a spring connecting two interfaces. Then Eq. (3) is no longer symmetric and 
Eq. (4) becomes 


<7i+l 

Qi 


Mi,6i- 1) 

1 


qi 

0 J U»-i. 



* = 1 JV, 


( 6 ) 


where S x and 3 x —i are random variables which correspond to the mistimed parameter on each side of interface 
z. Again, the frequency dependence of a and (3 is implied. Note the use of a semicolon in the index. 


2.1. Example: A two-degree of freedom per bay model. 



bay is shown with its two coupled neighboring blades. 


per 
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Let us illustrate the above with an example. Consider the mono-coupled blade assembly with two degrees 
of freedom per bay shown in Fig. 2. Dye and Henry [3] were among the first to propose this blade assembly 
model and it has been used subsequently both for blade assemblies [1] and for large space reflectors [15]. 
system will be utilized for demonstrative purposes for the remainder of the paper. In Fig. 2, represents the 

single-mode motion of the blade and qf accounts for the motion of the disk at the blade root. 
to q b are the blade (modal) mass m b and (modal) stiffness k' b . The mass m d simulates the effechve mass 
of the blade root and the corresponding section of the disk. The stiffness ^represents the basics 
of the rotor disk, whereas k' c provides disk coupling between neighboring blades. It is assum^that all 
bays have identical masses m b and m d and that the stiffnesses k' d , and k c may er m o y 
the next This reduces complexity and provides adequate means of mistunmg the natural frequencies of the 
bays. The average values of the stiffnesses are 4, k d , and k c , respectively, and correspond to the 
stiffnesses in a tuned assembly. We aim to study how different combinatmns of stiffness ^ m^sv^ires 
affect the sensitivity of the system to mistuning and how the sensitivity differs when the blade stiffness, the 
disk stiffness, or the coupling spring stiffness are mistimed. 

We now proceed to formulate a transfer matrix representation 
of the system in Fig. 2. Bay i is defined as the t* blade-disk element, 
including the spring connecting it to blade (i + 1) (see Fig. 3). The 
i th interface is the point where bay i and bay (i - 1) meet All 
parameters of bay i except the spring stiffness appear solely in 
relation with the I th interface, whereas k' c connects interfaces i and 
i + 1. For that reason, a dynamic force equilibrium at interface i 

will include both k[ and k\~ 1 . This yields the equations of motion | k‘ 



K 


rribq* + k&(<Zi — 9?) - 0, 


(7) 


^fcVVVV^ 


K- l (qf-qU)+K(qf -qi^qf+kHqf-q^+maqf ^0, ( 8 ) 

where Eq. (7) is local to the i* blade. We write the mistimed 
stiffnesses as 

fc* = k d ( 1 + Sf), k' b = k b ( 1 + ff) and fc‘ = k c ( 1 + 6f), (9) 


k' 



Figure 3 The i* bay. 


where S d 6 C and 6 b t are random variables with zero mean and standard deviation s d , s c and s b , respectively. 
Assuming harmonic motion, q b = -oj 2 q b , Eq. (7) may be rearranged as 


where 


rrihOJ 
u = h 


Q2 \ q b 

( 10 ) 



( 11 ) 


is a dimensionless frequency. Equation (10) is now used to eliminate the blade coordinates from Eq. (8), 
which becomes, 

-=..-,21 


d _ l 
" 1 + 6 ? 


k d 


tj 2 (l + S'? ) mil! 2 | d 1 + 6 f_ 1 d 


2 + S i +5 i-' + T c il+6i) ~ k c (l + %~* 2 ) ' k c J 9i 1 + *« C 




( 12 ) 


where we have introduced the following dimensionless parameters 

r k d l - kc ™ ^ 

k b k b irij> 


(13 a,b,c) 
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Defining 


and 


we can write 


pi;i - 1 = 


l+6f 


2 + 6; + 6f_ t + y^(l + $f) — j 


w 




- -2" 
mu 


k c (l + 6l-u> 2 ) k c 




j ±*U 

1 + 6 ? 


Qii 

Qi 


i 1 o J K-J’ ’ • 


(14) 


(15) 


(16) 


where the state vector for the system in Fig. 2 consists of the disk generalized coordinates. Equation (16) 
corresponds to Eqs. (5) and (6). However, we note that only in the case of coupling spring mistuning, i. e. , 
6f J 0, is the more complex form of Eq. (6) required, as otherwise aj ;i _i = 1 and = /?, . Thus, we have 
effectively formulated the dynamics of the system in terms of one reference degree of freedom per bay, in 
this case the coupling coordinate qf. Since the system is only mono-coupled, the other coordinate (the blade 
degree of freedom) has been eliminated. Of course, once the state vector has been calculated via transfer 
matrix methods, the blade dynamics are recovered from Eq. (10). 

In general, a mono-coupled assembly with P degrees of freedom per bay has P — 1 equations local to 
the bay, such as Eq. (10), and a single equation relating coupling coordinates of adjacent bays, e.g., Eq. (12). 
Needless to say, the complexity of this latter equation increases considerably as the number P is increased. 

For a tuned assembly we have 6? = 6f — = 0 and T 0 has the form in Eq. (4), where f3 a is given by 


Po 



-2 2 

uj mu* 

k c ( 1 -& 2 ) k c 


(17) 


In their work Cornwell and Bendiksen [15] and Bendiksen [1] took a different approach to the analysis of 
the same model. Using the theory of circulant matrices, they chose to eliminate the disk degrees of freedom 
by applying an approximate reduction procedure, arguing that in many cases of practical interest the disk 
stiffness is much greater than the blade stiffness. This is equivalent to setting = 0 and thus m = 0. 
We have found this to be an unnecessary approximation, although the approach allows one to evidence the 
coupling among all blades (as opposed to nearest neighbor coupling) through the disk degrees of freedom. 
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3. Propagation of Waves in the Tuned Structure. 

We examine the propagation of waves in a general, mono-coupled, periodic structure whose dynamics are 
governed by Eq. ( 4 ). It has long been known that energy-carrying motions in periodic structures only occur 
in isolated frequency ranges known as “passbands.” Outside the passbands only attenuated standing waves 
(or, more rarely, complex waves) can take place [16]. A physical understanding of these wave-propagation 
chara cteristics can be gained through the diagonalization of Eq. (4). This requires the solution of the 
eigenvalue problem: 

<Zi+ 1 
. ?< 


= T 0 


qi = A qi , i = 


(18) 


Equation (18) yields the eigenvalues and eigenvectors of T 0 . Its significance is as follows. The eigenvectors 
of T 0 define wave-modes, or characteristic waves, which propagate along the structure in such a way that the 
state-vector is multiplied by a complex scalar, A, as the wave passes through each bay. Thus the eigenvalues, 


a,j = T <»> 

define the frequency-dependent propagation properties of the corresponding wave-modes. The transfer 
matrix T 0 in Eq. ( 4 ) is real with detr o = 1, resulting in eigenvalues that are reciprocal and are either real 
or complex conjugates. We choose the convention that Ai has modulus greater than or equal to one. This 
associates Ai with the wave-mode traveling or attenuating in the direction of decreasing bay number. We 
shall call this direction the left or counterclockwise direction. The wave-modes appear as the eigenvectors, 
[1 ,A 2 ] t and [1, Ai] T , corresponding to the eigenvalues Ai and A 2 , respectively. The condition AjA 2 = 1 
indica tes that the two wave-modes could also be written as [Aj, 1] T and [1, Ai] T , which shows that they 
are equivalent except for their direction of travel. This is supported by the symmetry of the problem to 
clockwise or counterclockwise numbering of the bays. The independent wave-modes define the preferred 
means of wave propagation along the periodic assembly — much like normal modes are the preferred form 
of vibration of a structure. The normal modes form a basis for all vibration shapes in a structure. Similarly 
all possible waveforms in a structure may be written as a linear combination of the pair of wave-modes. 

The eigenvectors of T 0 are arranged as the columns of the matrix X : 

X = 

The matrix X defines the transformation 

Qi 

qi - 1 



from physical coordinates to left- and right- traveling wave coordinates at bay i, corresponding to the wave- 
mode basis. The displacement transfer matrix T 0 thus is transformed into a diagonal wave transfer matrix 


as 



L 

R . 

J l 


\r—Jrp Y T Aj 

- A 1 0 A |^J . - 



(22) 


Let Aj = define the complex propagation constant , p, with p = 7 + jo. Here 7 , the real part of the 
propagation constant, is the rate of exponential attenuation of the wave amplitude from one bay to the 
next. The imaginary part, a, is the interblade phase angle, the difference in phase between the motion of 
adjacent bays. In the literature dealing with wave propagation in periodic systems [16] the interblade phase 
angle is usually called the wave number. The propagation constant contains all the information about the 
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frequency-dependent propagation of waves through the assembly. Since /?„ in Eq. (19) is a real valued 
function of dimensionless frequency, u>, we distinguish between the follow ing cases: 

|A>(<*>)| < 2 In this case Ai and A 2 are complex conjugates of magnitude 1 , yielding 7 = 0 . These 
frequencies define a pass band in which waves travel without attenuation. From the real part 
of Eq. (19), the interblade phase angle is related to frequency by the dispersion relation: 

2 cos o = p o (u>), 0 < a < 7T. (23) 

For a given value of a, Eq. (23) has as many frequency solutions as there are degrees of 
freedom in each bay [16]. Hence, the number of passbands equals the number of degrees 
of freedom in each bay. One can reason that the natural frequencies of the system will 
lie in the passband because a perfectly periodic system must have periodic modes. In the 
blade assembly literature [1-9] these modes are usually referred to as constant interblade 
phase angle modes, i.e., modes with identical phase change from one bay to the next. The 
frequencies at which o — 7 r /2 will be referred to as midband frequencies. Note that the 
midband frequency is not necessarily located close to the mean frequency in the passband. 

\Poiu)\ > 2 The eigenvalues are real, hence 7 is nonzero, leading to attenuation. Also, adjacent bays 
are vibrating either in phase or out of phase, a = 0 or a = v, which implies that these are 
standing waves. These frequency ranges define stopbands. 

|/? 0 (£))| — 2 This gives 7 = 0 and o = 0 or o — tr, and defines the bounding frequencies or the passband/- 
stopband edges. 

3.1. Example. 



Figure 4 The passband/stopband structure of the assembly in Fig. 2, for k c — kj — ifi — 1 (tuned case). 

Note the two distinct passbands where 7 = 0 . Also note a singularity at ui = 1, corresponding 
to the natural frequency of a cantilevered blade. 
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Figure 4 depicts the passband/stopband structure for the assembly in Fig. 2 for one set of parameter values. 
The rate of exponential attenuation, 7 , and the inteiblade phase angle, a, are displayed as a function of the 
dimensionless frequency. As predicted above, two frequency passbands are observed, corresponding to the 
two degrees of freedom per bay. In the passbands 7=0 and unattenuated propagation of waves occurs. 
The wave travel is evidenced by the change of phase, 0 < a < 7 r, from bay to bay. The other regions are 
stopbands, 7 ^ 0, where standing waves decay exponentially. In stopbands, neighboring bays are either 
vibrating in phase or out of phase, hence there is no propagation. A special feature for this structure is the 
infini te attenuation observed at (D = 1, that is, at w 2 = ^. We explain this behavior by pointing out that 

this is the natural frequency of a blade cantilevered at its root, hence at this frequency qf = 0. Equation 
(10) confirms this. If the disk deflection equals zero, energy cannot be transmitted along the bays, hence the 
infinite attenuation of the waves, which can be viewed as an antiresonance. 


4. Natural Frequencies and Mode Shapes of Finite Cyclic Systems. 

We now make use of transfer matrices to derive a closed-form equation for the natural frequencies and 
normal modes of a finite, perfectly cyclic system. We already know from the theory of circulant matrices 
[17] that all cyclic structures share the same set of global modes, also referred to as constant interblade phase 
angle modes. We also know that structures described by real symmetric and circulant matrices have real 
eigenvalues, almost all of which are double. This is verified here using transfer matrices. 


4.1. Natural Frequencies of a l\med System 

In searching for the natural frequencies of the system whose dynamics are described by Eq. (4), we recall 
that (3 0 is a function of frequency, u>. Using the cyclic nature of the N- bay system and Eq. (4), we obtain 


92 

9i 


9JV+2 1 = t 0 n f 92 ] • 

9tf+iJ L9lJ 


(24) 


Hence nontrivial solutions are obtained if and only if 

det(T o "-j)=0. (25) 

Applying the appropriate transformation to T a yields a Jordan canonical form J ■ This allows for the case 
Aj = A 2 = 1, which does not have a complete set of eigenvectors. Thus Eq. (25) transforms to 

det {J N — J) = 0 . (26) 


Since the matrix J is an upper triangular matrix with the eigenvalues \\ and A 2 (from Eq. (19)) on the 
diagonal, J N is also upper triangular with Aj v and on the diagonal. Hence Eq. (26) becomes 

(A* -l)(Af -1) = 0, (27) 


or 


A lt 2 = Vl = e j2 * (n -' )/N n=l,...,N. 

Thus, from Eq. (19), the natural frequencies of the finite cyclic system, £>„, satisfy 


PoiUn) 


Poi&n) 


— 1 = e 




n — 1 , . . . , N, 


or 


( 2ir(n — 1 )\ . Ar 

Po&n) = 2 cos ( — J n = 1, ... ,1V. 


(28) 


(29) 

(30) 
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For each n, the number of frequency solutions to Eq. (30) equals the number of degrees of freedom per bay, 
hence an N- bay assembly with P degrees-of-freedom per bay has NP natural frequencies. An important 
consequence of Eq. (30) is /? 0 (w n ) = /? 0 (<*>.w-n+ 2 )> i-e., all natural frequencies are double, except u>\ and, 
for a system with an even number of bays, u>^+ q/ 2 - Note that the numbering scheme presented here does 
not place the natural frequencies in an ascending order. The relationship between the natural frequencies is: 
a>i < = wjv < o >3 = cD/v_i < . . . < d>(^+D /2 = £(.tv+i)/ 2 +i» if N is odd. If N is even the relationship 

is: Wi < u> 2 = u>N < • . • < &n / 2 = u N / 2+2 < u N f 2 + 1 , with unique frequency values at both ends. This 
pattern is repeated in each passband. The passband edges are natural frequencies for an A -blade system 
in the following cases. The lower passband edge is always a natural frequency corresponding to n = 1, 
or a = 0. If and only if N is even is the upper edge a natural frequency, corresponding to n = N/2, or 
a = 7T. Frequencies corresponding to all the other values of n are double. Furthermore, Eq. (30) gives 
-2 < Po(&n) < 2, indicating that the natural frequencies lie in the system’s passbands (see Eq. (23)). Thus, 
the assembly has N natural frequencies, mostly double, in each of its passbands. 

At this point, let us digress to examine briefly what would be the effect of damping and that of 
aerodynamic coupling (limited to adjacent blades), hitherto excluded. Equation (4) would instead be 


9«+i i r a>(<*>) 
. ?.• J i 



A>M, a c (w) G C. (31) 


The effect of damping is to render f3 0 complex, whereas aerodynamic terms cause p a to be complex in 
addition to introducing the complex off-diagonal term a 0 . The significance of a 0 is that the system is no 
longer symmetric with respect to clockwise and counterclockwise numbering. This makes sense since wave 
propagation in the aerodynamic system depends on whether the direction of wave travel is the same as the 
direction of rotor rotation, or not. The eigenvalues of 1 0 are then 



Po,a a e C, 


(32) 


instead of the values expressed in Eq. (19), with AjA 2 = a 0 . The fact that (3 0 is complex precludes the 
existence of passbands, as expected in a nonconservative system. Equation (30) for the natural frequencies 
takes on the more complicated form 


&(£>») = e i 2 *< n -'V N +a 0 (u) n )e-^ ( "- , >/ JV 


n = 1 . . . N, 


(33) 


indicating that the presence of a 0 , which is due to aerodynamic terms, splits the natural frequencies of the 
system. 


4.2. Nodal Diameter and Nodal Circle Modes. 


From Eq. (30) there are N natural frequencies in each passband and we wish to determine the corresponding 
modes shapes. The complex mode shapes can be found simply through a recursive application of the wave 
modes given in Section 3. This yields the shape [1, A,, 2) A^, . . . , A*" 1 ], with A u given by Eq. (28). The 
mode shape may also be given in terms of the interblade phase angle for the mode. A comparison of Eqs. (23) 

and (30) shows that to each u„* corresponds a mode shape with constant interblade phase angle <r„ given 
by 


_ 27r(n — 1) 
“ N ’ 


n=l...N. 


(34) 


This result is independent of the passband number and implies that for a given n all frequency solutions of 
Eq. (30) are associated with the same global mode shape, characterized by a constant interblade phase angle 


* Recall that the index n on natural frequencies does not imply an ascending ordering of natural 
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along the assembly. By global we mean that modes of interblade phase angle a n , in the different passbands, 
feature the same deflection pattern from bay to bay along the assembly and differ only by their local vibration 
shape within an individual bay. In the first passband all blade-disk elements vibrate in their first mode, in 
the second passband in their second mode, etc. From Eq. (34) we note that expO'a n ) = exp(- jff(jv+2-n)X 
which means that to each double natural frequency, u>„ = d>(jv+2-n)» corresponds a complex conjugate mode 
pair, one traveling clockwise and the other counterclockwise. Since they occur at the same frequency, these 
two traveling mode shapes may be combined to form two orthogonal real (standing waves) modes, simply 
by separating their real and imaginary parts. 

The global modes are often classified in terms of their interblade phase angle. The first mode would 
thus be a zero interblade phase angle mode, with all bays moving in phase. The second mode and mode N 
are interblade phase angle modes, etc. Another means of classification is to consider the number of nodal 
diameters featured by the mode. A nodal diameter is a term borrowed from continuous two-dimensional 
systems with cyclic symmetry, where it is a nodal line passing through the center of symmetry. In our case 
it is simply a diameter connecting pairs of nodes or bays with zero deflection positioned on opposite sides 
of the structure. The first mode is thus a zero nodal diameter mode. The second mode and mode N are 
one nodal diameter modes, etc. Another concept from continuous two-dimensional cyclic systems is that 
of no dal circles. Here, a nodal circle is drawn through the nodes of the individual bays of the assembly. 
Naturally in a cyclic system all N blades have nodes at the same distance from the center of the rotor. We use 
the number of these nodal circles to distinguish between modes with the same global shape but in different 
passbands. In general, in the first passband the blade-disk elements have no nodes, hence all modes in the 
first passband are zero nodal circle modes. In the second passband all blade-disk elements typically feature 
one node, hence all second passband modes are one nodal circle modes, and so on. 



Figure 5 The global mode shapes of a 28-bay cyclic assembly. The 1 st mode is a unique zero nodal 
diameter (o\ = 0) mode. The 2 nd mode is a one nodal diameter (cr 2 = fj) mode. The 28 
mode (not shown) also has one nodal diameter but with <728 = -§• The 15 th mode is a 
fourteen nodal diameter (<715 = 7r) mode and is unique. 


Figure 5 illustrates the global modes of a 28-bay cyclic system. For modes corresponding to the same 
double natural frequency only one is shown. For example the 28 Ul global mode would be like the 2 global 
mode except for the direction of travel. Combining global modes 2 and 28 would yield two standing mode 
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shapes, identical to the one shown with one mode rotated 90 degrees with respect to the other. 

4.3. Example. 


7 
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<3 2 

3 


2 

1 


0 

0 N/4 N/2 

number of nodal diameters 

Figure 6 Natural frequencies as a function of the number of nodal diameters, for both passbands of the 
system in Fig. 2 with k c = k d = rh = 1. 



Turning our attention to the system of Fig. 2, we have, from Eqs. (30) and (17) 


k d uil mu 

P(&n) = 2+J±- - --T 

&c( 1 k 


mu>n 

^-2- = 2 cos 


27r(n — 1) 
N 


which when solved for yields 



where the minus sign yields values in the first passband and the plus sign gives second passband frequencies. 
In Eq. (36), a n is the interblade phase angle of the n lh mode in the passband, given in Eq. (34). 

It is customary to plot the natural frequency distribution against the number of modal diameters. This 
gives an indication of the number of natural frequencies per unit frequency, the modal density, in a given 
frequency range. In Figure 6 this is done for m = 1, k d = 1 and k c = 1. An interesting feature, sometimes 
overlooked, is that the curves have local extrema at zero and N/2, hence the modal density is infinite at 
passband edges. The curves increase monotonically from 0 to N/2 and possess no other extrema. 
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I st passband 2 nd passband 

Figure 7 Typical local mode shapes in each passband with zero and one nodal circle, 


respectively. 


While all circulant systems feature the same global mode shapes, illustrated in Fig. 5, the shape local 
to each bay is determined from Eq. (10), which may be used to calculate the radius of the nodal circle for 
the mode. For this system the modes of the first passband (d; 2 < 1) have no nodal circle (q d and q b have 
the same sign) and the modes of the second passband (<D 2 > 1) have one nodal circle (q d and q b have the 
opposite sign). From Eq. (10) we also see that the ratio between q d and q b is frequency dependent, hence the 
nodal radius varies within the passband. The local vibration patterns of the blade-disk elements are shown 
in Fig. 7 for the two passbands. 

Appendix A contains a detailed discussion of the effect of the system parameters on the width and 
location of the two passbands. The main findings are that the lower passband edges correspond to a mode 
which has all blades vibrating in phase and thus to a system in which the blades are completely uncoupled. 
By that we mean that no interaction occurs between bays. Also, a system is called weakly coupled in some 
passband if the passband is narrow. 

5. Natural Modes of a Mistimed Structure. 

Let us begin the analysis of mistuned systems by examining the effect mistiming has on the natural frequencies 
and mode shapes. When mistuning is introduced, the system is no longer symmetric with respect to the 
clockwise and counterclockwise directions. This loss of cyclicity results in a splitting of the double natural 
frequencies, such that the system has P clusters of N distinct natural frequencies, where N is the number 
of blades and P is the number of degrees of freedom per bay. The clusters of frequencies correspond, 
approximately, to the passbands of the system’s tuned counterpart, although they generally are wider. The 
corresponding mode shapes are standing waves that no longer possess the cyclic symmetry exhibited by 
the tuned system, where all blades vibrate with the same amplitude. Instead, the vibration energy may be 
concentrated in a handful of blades that have significantly larger deflection than the majority of blades. We 
refer to this phenomenon as mode localization. 

As an example of the effect of random mistuning, we consider the mode shapes of the mistuned blade 
assembly model depicted in Fig. 2. Figure 8 displays the global blade-to-blade pattern of the mode shapes 
in the second frequency cluster of a 28-blade assembly, where the individual blades vibrate primarily in 
their second mode (see Fig. 7). High levels of localization are observed. The mode shapes have changed 
drastically compared to the tuned modes depicted in Fig. 5, and the vibration energy is no longer uniformly 
distributed along the rotor, but concentrated in as few as 5 blade-disk elements. This may result in a 
substantial increase in amplitudes and a potential decrease in fatigue life. 

A compact characterization of the strength of localization and the prediction of high sensitivity to 
mistuning is the subject of the remainder of this study. 
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Figure 8 Localized modes in the second frequency cluster of a 28 blade assembly with k d = m = 1, 
k c =0.1, and uniform random mistuning in the blade stiffness with 5% standard deviation 
(sb = 0.05). Modes are sorted by increasing natural frequency. 
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Figure 9 Scattering of waves at substructure interfaces. 


6. Waves in Mistuned Systems. 

A statistical approach is chosen to examine the high sensitivity to mistuning and the occurrence of localization 
in randomly mistuned assemblies. A statistical description of the modes of vibration is not practical due to 
the switching of the associated natural frequencies as the strength or the distribution of mistuning varies. 
We choose instead to examine the propagation of incident waves in mistuned assemblies. This allows us 
to control the frequency at which we wish to examine localization. In order to avoid the contamination of 
the localization effect by the cyclicity condition, an infinite assembly is studied. We further assume that 
mistuning is restricted to a segment of N bays (numbered from i = 1 to i = N) embedded in an otherwise 
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tuned infinite assembly. We are interested in the transmission of incident waves along the mistimed segment. 
An advantage of this approach is that a traveling wave that exits the mistimed segment will be propagated 
away and will not return. 

In Section 3 we demonstrated how the transfer matrix T 0 of a tuned bay is diagonalized by the 
similarity transformation defined by the matrix X , which has as its columns the eigenvectors of T 0 . We 
called the diagonalized matrix a wave transfer matrix, W a . The transformation introduces a new set of 
coordinates, namely the left- and right-traveling components of a wave. The transformation defined by X 
does not, in general, diagonalize the transfer matrix of a bay belonging to the mistuned segment, T { or 
(i = 1, . . . , N), as the case may be. Instead a non-diagonal wave transfer matrix is generated. The elements 
of this matrix may be expressed in terms of the more familiar reflection and transmission coefficients, as 
follows. 

At the interface between dissimilar bays, waves are split into a transmitted part and a reflected part. 
As illustrated in Fig. 9, the left traveling wave incident to bay i — 1, i,, is the sum of a transmitted left 
traveling wave, t,L ,+ 1 and a reflected right traveling wave, riRi. Likewise, the right-traveling wave incident 
to bay i + 1, ii,+i, is comprised of a transmitted part, t, R , , and a reflected part, . Here <,• and t, are 

transmission coefficients and r, and f; are reflection coefficients in the left and right directions, respectively. 
Transmission and reflection coefficients are the complex amplitudes of transmitted and reflected waves due 
to an incident wave of unit amplitude, respectively. The above defines a scattering matrix. Si. 


Li 



(37) 


where directional symmetry in the absence of aerodynamic forces dictates that U - U and f; - r,-. Solving 
Eq. (37) for R i+ 1 and L i+ \ yields 


r 1 -Ti 1 



Equation (38) tells us that the off-diagonal elements in W t govern what portion of a right or left traveling 
wave is reflected, thereby generating a wave in the opposite direction. The transmitted portion of the incident 
wave is determined from the first diagonal element in the wave transfer matrix. For a tuned bay W, = W D 
is diagonal and there is no reflection at the interface. Hence r< = 0 and the wave is fully transmitted. For a 
mistuned bay there is a reflection, or scattering at the interfaces between bays. A wave incident to a segment 
of randomly mistuned bays will experience multiple reflections whose effect may be to trap a wave near the 
incidence region. Only a frequency-dependent fraction of an incident wave is transmitted along to the far 
end of the mistuned segment. This effect is called localization. In a mistuned system passbands no longer 
exist since all waves are attenuated. 

For a segment of N dissimilar bays the wave transfer matrix is the product of the random wave transfer 
matrices of the individual bays 


jv 


r 1 


w N = n = 


~PN 

tn 


t n t n ^ | (39) 

PN ~ 

. Tn ~ i 

L Tjf T N J 

The one-one term, l/r^r, tells us which portion of an incident wave is transmitted to the far end of the 
mistuned segment. Obviously the behavior of the transmission coefficient for the mistuned segment, r jV , 
governs the strength of the effects of mistuning and the resulting localization. We define the localization 

factor as 


7 = lim 

N->a 


1 . 1 1 

77 In — , 

' L N Tjy J 


(40) 
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implying that asymptotically, the ratio of emergent to incident wave decreases exponentially with an increas- 
ing number of bays, N . The localization factor, 7 , defines the average exponential decay rate per bay and 
thus is a descriptor of the strength of localization. Assuming that the random variables <5, (i = 1 ... N) from 
Eqs. (5) and ( 6 ) form an ergodic sequence, then the transmission coefficient, r^(5j, . ..,6 N ), is alsoergodic. 
Ergodicity implies that, with probability 1, the limit in Eq. (40) is equivalent to the ensemble average for a 
finite mistuned segment of length N: 


/ 1 

1 

( 


\ 2 N 

tn 


where () denotes the expected value of a random variable. 

The assumption that a cyclic structure is infinite is of course a limitation. It must be understood that this 
type of analysis is only applicable if the localization is strong enough that the vibration energy is confined 
to a sufficiently small region compared to the size of the assembly. For a nearly cyclic finite system this 
basically requires that the wave be adequately attenuated before reaching the incidence region again and 
thus interacting with other waves from the same source. 

6.1. Sensitivity to Mistiming. 

In general Eqs. (40) and (41) for the localization factor cannot be evaluated in closed form. However, 
analytical approximations are possible using perturbation methods, as follows. 

We seek an expansion of ^ in Eq. (40) in terms of the small mistiming parameter, 6 ,- (8, of order e or 
smaller), thus treating the mistuned system as a perturbation of the tuned system. This is obtained through 
a Taylor expansion of the function /?(<$,) in Eq. (5): 

m) = f3 0 +p'(p)6 i+ 0(6}) (42) 

and subsequent expansion of t n and expansion and averaging of 7 . The expansion in Eq. (42) must of 
course be uniform, i.e., (3'(0) must be of order one (order e®) or smaller. However, if the expansion in 
Eq. (42) is nonuniform, i.e., when is large (order e -1 or larger) the technique breaks down. In this case 
the perturbation expansion is only valid for 6, second order (order e 2 ) or smaller. This breakdown indicates 
high sensitivity to mistuning and it is when the breakdown occurs (if it occurs) that systems have been seen 
to enter the realm of strong localization. This suggests the use of the first-order Taylor coefficient of /?(£, ) 
as a measure of sensitivity to mistuning: 

S = (3'( 0 ) (43) 

A similar scenario may take place in the Taylor expansion of the transfer matrix in Eq. ( 6 ), albeit with four 
different first-order Taylor coefficients for a and /? with respect to <5, and £_i. We do not formally define a 
sensitivity measure for this case. 

When the sensitivity measure, S = 0’( 0), becomes large the expansion in Eq. (42) is nonuniform, 
indicating a qualitative change in the assembly’s dynamics, and other avenues must be explored. Note that 
in the passbands of the tuned system we have -2 < (3 0 < 2, but when S is large, (3{8 t ) in Eq. (42) has 
the potential to become large for first-order mistuning in that frequency range. This suggests treating the 
off-diagonal terms in the transfer matrix of Eq. (5) as small perturbations compared to the large /3(8 t ), thus 
yielding an alternative expansion for the high sensitivity case. 

The two techniques outlined here, for the cases of normal and high sensitivity, are detailed in the 
following sections. 
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6.2. Classical Perturbation Method — Normal Sensitivity. 

Consider the case where the expansion (42) is valid for first-order mistiming (ft of order e), i.e., when 
5 = ft(0) is not large ( S of order one or smaller). We shall refer to this case as low or normal sensitivity. The 
much simpler case in Eq. (5) is treated first to better explain the technique. In Eq. (5), ft = ftft) depends 
only on one mistuned parameter, with /?(0) = ft corresponding to the tuned case. A uniform expansion of 
Ti in the small mistuning parameter ft yields 


Ti = 


m) -i 
1 0 


ft - 

1 0 


i"L r/m oi, r/m oi g 
I J [ 0 OJ [ 0 OJ 2 "■ 


. 6 } 


(44) 


= To + T'ft + T'y + . . . 


Ot 


(45) 


A tr ansf ormation to wave coordinates is accomplished with the matrix of eigenvectors of T, 

X - ^ ! 

[ e ~ ]<T e ltr 

where c is defined in Eq. (23). Note that the analysis is limited to the tuned system’s passband, since there 
is already large attenuation in the stopband. The wave transfer matrix is expanded as 

(46) 


Wi = X~ J TiX = W 0 + W'ft + W" y + • ■ 
t* 9 0 1 [ft* 1 1 ftft( 0) 

0 e -J<T + -1 e~ j<T 2 ? sin a 


ft' 1 ] <S 2 ft'( 0 ) 

[— 1 e~i a J 4 j sincr 


We need to evaluate the wave transfer matrix for a segment of N mistuned bays, each with its own random 
mistuning parameter ft. The random variables ft corresponding to each bay are considered independent 
and identically distributed, with zero mean and standard deviation s*. An assembly of N bays has the wave 
transfer matrix Wjv, which, when expanded to the second order in the ft’s, becomes 


- n ( w ° + s ‘ w * f w ") + °tf> 

i=N ' ' 

=W 0 - v +f; Wj-'W' + f; wj-' W" wg- > §.££ R Mi + 0(6]). 


(47) 


The matrix R in Eq. (47) has a highly complex form, not evaluated since that entire term will vanish in 
the averaging process that follows, due to independence of the random variables ft, (i = For an 

approximation of the localization factor, only the first diagonal element of War, tuff = 7 * 7 . is needed. We 
have 


J- = ft*' 

tn 


1 + 


ft(0) 


N 




$"(0) 


N 




2 j sin a * 4 j sin a 

i=l t=l 


- 0 ( 6 ?), 


1 

2 

TN 



N 


N 




+ 0 ( 6 |). 


(48) 


(49) 


4sin2(7 ,=i i=i J 

Since the ft’s have zero mean and are independent, we have (ft ft) = 0 for j ^ i and (ftft) = s 2 . Using 
the expansion log(l + x) = x + 0 (x 2 ), we obtain the approximation of the localization factor, valid for low 
sensitivity, 

2 


7 = 


2JV 


In 


tjv 


^ 1 / ft (0)5 A 2 = (5 (o>)) 2 5 2 

— 2 \ 2 sin a / 2(2 + ft(w)) (2 — 


ft(a>))’ 


0(S ) < 1 , 


(50) 
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where S is the sensitivity measure defined in Eq. (43). Note that both S and /3 0 are functions of u> and vary 
within the passband. The localization factor, 7, allows for localization effects to be characterized in a simple 
and compact way, without extensive simulations of mistuned systems. Note that Eq. (50) is general and 
applies to any mono-coupled mistuned assembly, as long as the mistuned parameter is not one that connects 
the two interfaces bounding the bay. It is of interest to note that the second-order Taylor coefficient, /3"(0), 
is cancelled from the derivation. According to Eq. (50), the onset of localization increases as the sensitivity 
measure squared. Also, to first order, the localization factor increases with the square of the mistuning 
standard deviation. The approximation of the localization factor becomes unbounded at frequencies which 
correspond to the edges of the passband of the tuned structure, 0 O = ±2. This is reasonable, since at those 
frequencies the transformation matrix X (see Eq. (45)) is singular. 

In order to obtain a quick estimate of the level of localization present in a system, it is advantageous 
to focus on the localization factor at the midband frequency. This frequency corresponds to the interblade 
phase angle o = 7t/2, thus is the median natural frequency* of the tuned system. As seen in Appendix A, the 
midband frequency need not be located close to the mean frequency of the passband. One has, at midband. 


* Ymid — 




8 


(51) 


which gives a good indication of the strength of localization effects in a “typical” mode of the system. 


Let us now consider the case where a parameter common to two adjacent bays is disordered. In this 
case the transfer matrix has the general form in Eq. (6), where a and (3 are functions of two independent and 
identically distributed random variables, 8 t and $,_i. Assuming that /? and a may be uniformly expanded 
in the small disorder parameters <5, and <5,_i, we find an expansion of 
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=T 0 + T,jSi + T, 2 Si - 1 + + T,i 2 SiSj + T, 22 -^- + ..., 
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0 

+ 


(52) 


where the index notation T \ and T f 2 has been adopted to denote differentiation with respect to the first variable 
(Si) and the second variable (6 t _i), respectively. A transformation to wave coordinates is accomplished 
using the matrix X from Eq. (45). The result is an expansion of the wave transfer matrix in the small disorder 
parameters Si and _ 1 : 


= W 0 + W.,6, + *W, U 6 1 + W, l2 t,s,^ * w, 22 ^- + ... 

The wave transfer matrix of an A -bay segment is 

Wiv = W N . N _iW N - hN-2 ••• ^2;J% 

- II ( + + W ’2*i - 1 + + W, 12 SiSi. , + W, 22 ^ . 


(53) 


(54) 


<7 = 7r/2 does not correspond to a natural frequency unless N/4 is an integer. Nonetheless, in each passband, the frequency corresponding to 
<7 - tt/ 2 (see Eq. (23)) will have an equal number of frequencies above and below it 
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Carrying the product out to the second order in the 6’ s, the one-one term is found to be, after tedious algebra. 
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(55) 


, Ci,C 2 €C 


The tilde on underlines the fact that a few liberties have been taken, i.e„ we have ignored very complicated 

terms that we know will vanish in the averaging process at a later stage in the derivation. The term containing 
W ,2 W,i appears due to the occurrence of the same random variable in two adjacent matrices. Next the 
magnitude of is evaluated: 
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i 2 = 1 + iii 6i + Ri 22 6 < + R i’ R * elR - 


(56) 


1=1 


i=l 


Due to the terms linear in 6, a second-order expansion of log|ti>ff| 2 is required, that is, log(l + x)~ 
x - x 2 /2 + 0(x 3 ). Dividing by 2 N and taking the limit as N -* 00 yields the localization factor, after 

substantial algebra: 
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(57) 


(58) 


which reduces to the form of Eq. (50) for /?(£,, Si- 1 ) = /?($,) and a(6 t , 5,_i) = 1. Finally, in the numerous 
cases where f(f . . 

a(5,^_i) = ifei (59) 


m) ’ 


Eq. (57) simplifies to 


7 = 4 




(60) 


At midband Eq. (60) reduces even further, to: 
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6.3. Modified Perturbation Method — High Sensitivity. 

As suggested in Section 6. 1 and observed in Section 6 . 2 , the Classical Perturbation Method fails in the limit 
of high sensitivity, i.e., when 5 = /?'( 0) is large. Not only may this be seen by the failure of the Taylor 
expansion of 0 (Eq. (42)), but also by the fact that 7 in Eq. (50) becomes large (order e -1 ) for 5 large, while 
it supposedly is a small perturbation of its zero value in the tuned system’s passbands. The introduction 
of first-order mistuning in a system with strong sensitivity (/?'( 0 ) » 1 ) has the effect that in the passbands 
P(6i) becomes large 1 . This suggests the following modified perturbation scheme, where the off-diagonal 
terms are considered as the perturbation. 


T 1 . _ P($i) — 1 
1 1 0 


m ) 01 ro -11 
0 oj [1 oj 


= Tj + AT. 


(62) 


A brief review of the equation of motion, Eq. (3), aids in a physical interpretation of the Modified Perturbation 
Method. The larger 0 becomes compared to 1 , the less q { is influenced by q M and Consequently, 
in Eq. (62), a coupled system is treated as a perturbation of the uncoupled system. Since the modified 
unperturbed matrix T, is already diagonal, there is no need for a wave coordinate transformation. Hence 
the physical coordinates are also “ wave coordinates 2 " and J, is a wave transfer matrix for the unperturbed 
structure. The wave transfer matrix for the perturbed, N - bay structure is then 


w'-nu-nn+f: 

i=N i=N /=! 


where the product is expanded to the first order in the perturbation, AT. Due to the special form of AT all 
terms in the sum vanish except the first one and the last one. These terms become the off-diagonal terms 
in the matrix product and do not affect the transmission coefficient. Thus, in accordance with Eq. (41) the 
localization factor becomes 

7 = (In |/?(<5)|) = J In |/3(<5)|pdf 6 (,5)d<5, ( 64 ) 

where () denotes an average and pdf 6 (<5) is the probability density function of the random variable 6. 

Next we attempt to express the localization factor in terms of the large sensitivity measure, S = 0'(O). 
Focusing on the midband frequency, where the intetblade phase angle is o = tt/ 2 and (3(0) = 0, the first-order 
approximation of 0 mid for small 6 is (assuming that the remainder of the expansion is uniform) 
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(65) 


Equation (64) may now be written as 


7 mid - In |S TOi<1 | + J In |^|pdf { ( 6 )d£, ( 66 ) 

which in the case of a uniform distribution, 6 £ [~W, W] = [— \/3ss, \/ 3 s^ ], reduces further, to 

Imid ^ In |5mid| + ln(v^Si) - 1, 0(S mid ) > 1, (67) 


1 Recall that in a tuned system’s passband —2 < 0 ( 0 ) < 2. 

Since the unperturbed structure is uncoupled it does not actually transmit waves. 
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a remarkably simple and general characterization of localization in the high sensitivity case. Note that only 
the system sensitivity and the mistiming standard deviation are needed to evaluate 7 . Alsonote that since 
S a is large 7 w is not first order but of order one, corresponding to a strong localization behavior. Pierre 
[14] examined a chain of single-degree of freedom oscillators, where S = l/R, i.e., the sensitivity is the 
inverse of the dimensionless coupling between oscillators. His result verifies Eq. (67). 

In the limit of strong sensitivity, the Modified Perturbation Method treats a coupled system as a 
perturbation of an uncoupled, disordered system. This implies that the method is applicable only to weakly 
coupled assemblies and leads to an interpretation of sensitivity as an inverse of the coupling among bays. High 
sensitivity to mistuning has already been observed in systems with weakly coupled bays [14]. However, the 
definition of coupling in complicated multi-parameter systems is not always possible, whereas the evaluation 

of S = 0' (0) is straightforward. . . 

Note that the results of this section are valid only in the cases where the mistuning is not common 
adjacent bays as in Eq. ( 6 ). Since the systems of immediate interest to us do not exhibit strong localization 
when common parameters are disordered, this extension is not developed here. 



Figure 10 Approximation of the localization factor at the midband frequency as a function of the 

sensitivity measure (S mid «#,*(<»). ( ) corresponds to the classical perturbation 

result (Eq. (50)), valid for S mid small. ( ) corresponds to the modified approach 

(Eq. (67)), valid for S m id large. The mistuning is uniform with 5% standard deviation. 

Figure 10 illustrates the perturbation results for the mid-passband localization factor (Eqs. (51) and 
(67)) for a generic system described by Eq. (5), plotted as a function of sensitivity. The mistuning is uniform 
with 5% standard deviation. Equations (51) and (67) are valid in the limit of weak and strong sensitivity, 
respectively. Observe the rapid onset of localization, followed by a more moderate increase but large values 
of 7 mid in the limit of high sensitivity. 


22 




6.4. Example. 

We now apply the tools developed above to the system depicted in Fig 2. Corresponding to the mistimed 
parameters k' b and k' d are measures of sensitivity to mistuning, S blade and S disk , according to Eq. (43). No 
such measure has been defined when the mistuned parameter is one connecting two interfaces, therefore 
the case where k c is mistuned is given a different treatment. A parametric study identifies combinations of 
parameter values leading to high sensitivity. Analytical approximations of localization factors are obtained 
in the limits of weak and strong sensitivity. These approximations are verified by Monte Carlo simulations. 

6.4.1. Measures of Sensitivity 

Disk Mistuning: From the definition of 0 in Eq. (14), when 6f = 6? = 0, one has: 


k d , 




/*«?)« 2 + f(l+<f)- r _ 

kc K c (l — OJ 2 ) 


nujj 2 

IT' 


Thus, from Eq. (43), the measure of sensitivity to disk stiffness mistuning is simply. 


( 68 ) 


Note that since 0(6?) is linear with respect to the mistuning, its first-order Taylor expansion is exact. S disk 
is easy to comprehend. The sensitivity increases linearly with the disk stiffness and is inversely proportional 
to the stiffness of the coupling spring. It is independent of the blade stiffness and the mass of the blade 
and the disk. The sensitivity measure is also independent of frequency and is therefore the same in both 
passbands. Note the obvious correlation between weak interblade coupling and high sensitivity as k c — * oo 

and k d — 0. These results confirm those of Wei and Pierre [9], Pierre [14], and Cornwell and Bendiksen 
[ 15 ]. 

Blade Mistuning: For blade mistuning only, 6? = 6? =0 and Eq. (14) simplifies to 


0 ( 6 -) = 2 + ^ 1 _ TWfZ)2 
k c k c ( 1 + ^ — u> 2 ) k c 


(70) 


Thus, from Eq. (43) 


gblade = 

~ (1 — <jj 2 ) 2 k c 


(71) 


Equation (71) is more complicated than it first would seem, because the location of the passbands, and 
hence the corresponding ranges for a), depend on the system parameters k c , k d and m. The sensitivity is 
of course very different in the two passbands of the system. To give an example of how deceptive Eq. (71) 
is, one might wrongly assume that a large coupling spring stiffness, k c , always leads to low sensitivity. A 
counterexample is given below. However, small values of ic c do ensure high sensitivity, as we soon will see. 

Focusing on the midband frequencies, we solve Eq. (35) for w with o = 7t/2, and obtain 


+ k d + m + 1 y [2^ c + k d + m + l] 2 — 4m + 2Jt c l 
" 2r7i ’ ™ 

where the minus sign denotes the first passband. Substituting the values of into Eq. (71) yields 


2 _ 2 kc ^ 

^ mid * - i 
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- 4 Z 


2k c + k d - m - 1 q: 


\J [2k c + + m + l] 2 - 4m [k d + 2k c ] 


(73) 
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Contours of the sensitivity measure, , at the two passband centers. In each plot one of 
the parameters k c , k d and m is fixed while the others are varied. 














where the plus sign now corresponds to the first passband and the minus sign denotes second passband 
sensitivity. The sensitivity to blade mistuning at midband is plotted in Fig. 1 1. The plots show contours of 
the surface Sf£$ e (k c , k d , rh) by, in turn, fixing one parameter and varying the other two. A study of the 
plots yields information about the combination of parameter values that lead to strong sensitivity in the two 
passbands. 

Figures 11(c) and 11(e) tell us that strong sensitivity in the first passband may be produced by large 
as well as small values of k c , the coupling spring stiffness. At first it seems surprising that a large k c leads 
to high sensitivity at the first midband. Indeed, Fig. A2(c) in Appendix A indicates that for high values of 
k c the passband width is large. However, we see that for large k c all natural frequencies except the first 
one merge with the upper passband edge. Hence, although the passband is large, the range of the natural 
frequencies becomes small, and this system has very weak coupling. In Appendix A a system with stiff 
coupling springs is likened to a collection of blades mounted on the perimeter of a rigid ring. The blades 
are thus decoupled in the limit k c — ► oo, vibrating in their first mode with qf ~ 0. The opposite occurs in 
the second passband, where the bay is vibrating in its second mode that, unlike the first mode, requires a 
nonzero qf. In this case the stiff coupling spring contributes heavily to the coupling between bays and the 
sensitivity decreases as k c increases. 

In addition to the above observations we note from Fig. 1 1 that if m is large enough, increasing it 
will cause a decrease in first passband sensitivity and an increase in second passband sensitivity. Varying 
rh for small m has only a small effect on sensitivity, as evidenced by the lower portion of Fig. 11(e) and 
Fig. 1 1(f). Similarly we see that an increase in k d increases first passband sensitivity except when k d is small 
(Fig. 11(a)) or when k c is very large (Fig. 11(c)). Second passband sensitivity decreases with increased 
k d but only in ranges determined by m and k c . A combination of a small k c and a large k d leads to high 
sensitivity in both passbands. 

Coupling Spring Mistuning: Although we have not defined a sensitivity measure for the system in 
Eq. (6), we attempt to examine its sensitivity to mistuning by considering all four possible first-order Taylor 
coefficients. From Eq. (14), with 6f = Sf = 0, 
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(76) 

(77) 


All four Taylor coefficients in Eqs. (76) and Eqs. (77) remain of order one or smaller for all parameter values. 
We conclude that in the case of coupling spring mistuning the classical perturbation is valid for all parameter 
values. This system is thus never strongly sensitive to coupling spring mistuning and only weak localization 
may be expected to occur. 
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6.4.2. Localization Factors and Monte Carlo Simulation. 

Based on the expressions for the sensitivity measure, we are now in the position to examine the localization 
factor due to the above three types of disorder in the limits of low and high sensitivity. Using Monte Carlo 
simulations we confirm the validity of the perturbation approximations of the localization factor and the 
validity of 5 as a measure of sensitivity. 

Monte Carlo simulations were performed in the following manner. A series of random transfer matrices 
representing the bays of a finite, cyclic assembly were generated, based on a random sequence, <5,-, from a 
random number generator. We labeled the resulting assembly one realization of mistuning. The matrices 
were multiplied together and transformed into wave coordinates, as explained in Section 6, yielding the 
wave transfer matrix for the mistuned assembly (see Eq. (39)). The localization factor for this particular 
realization of an N -bay system may be calculated from. 




tn 


(78) 


This was repeated for a large number of realizations of the assembly and the result averaged, adding 
realizations until the desired accuracy had been reached. The above constitutes one Monte Carlo simulation, 
for a system with some set of parameters k c , kd and fh vibrating at some frequency. 

One point must be made about N, the number of bays in each realization. A product of a finite number 
of transfer matrices does not account for the infinite number of successive reflections that occur in the bays 
of an infinite system. A Monte Carlo simulation based on single bay realizations, N = 1 , would only account 
for direct transmission through the bay and all reflections would be truncated. A Monte Carlo simulation 
using two bay realizations, N = 2, would additionally account for the wave fractions that are reflected twice 
as they propagate through the assembly, but all higher order reflections would be truncated, and so on. The 
number of bays in each realization is especially important when localization is strong (the off diagonal terms 
in W, are no longer very small). For the cases studied, it was found that fewer than eight bays (N = 8) 
for each realization would, for strong localization, converge to an incorrect value for the localization factor. 
The number of realizations required for the convergence of 7 jv varied but was usually in the thousands. 

Verification by Monte Carlo simulations consists of two parts. On the one hand we verify the variation 
of 7 as a function of frequency by running a series of simulations for a range of frequencies and plot the 
results along with the perturbation approximations in the passbands. The other part of the simulation process 
is the verification of the two perturbation solutions of 7mtd as a function of the sensitivity measure, S m id- It 
is clear that both 7 mid and S mid are functions of k c , k d and m and may thus be parameterized with any or all 
of these parameters. This offers the opportunity to examine the invariance of the relationship between 7 m:< j 
and S m id with respect to the system parameters, such that we would expect to obtain the same relationship 
no matter which of k c , k d or m is varied. This relation should also be independent of which passband is 
chosen for the simulation. If this invariance is verified then the localization factor in Eqs. (51) and (67) will 
prove to be an extremely general tool for predicting mistuning effects when used in conjunction with the 
measure of sensitivity, S m ,d. 

For the Monte Carlo simulations a uniform distribution of width 2 W was assumed for the mistuning. 
The standard deviation of the mistuning is 


Disk Mistuning: 

From the perturbation analysis, the localization factors are: 
Weak Sensitivity 

s 2 d k\ 

<y, ~ - — 

8 fc 2 sin 2 o ’ 

s\k 2 d 

7d,w - -g fcT- 


(79) 
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Figure 12 Monte Carlo simulations ( ) are used to verify the invariance of Eq. (80) ( ) 

and Eq. (82) ( ) . The simulation result is a collection of data generated by varying 

both k c and k d in both passbands, with all results falling on the same curve. The standard 
deviation of uniform disk stiffness mistuning is 5%. 



Figure 13 Exponential decay, y d , due to 5% disk stiffness mistuning vs. frequency. In stopbands 

Monte Carlo simulations ( ) nearly coincide with the propagation constant for the 

tuned system ( ) . In (a) the simulations results overlap the classical perturbation 

result ( ) , Eq. (79), in the passbands. In (b) the simulation agrees with the modified 

perturbation result ( ) , Eq. (81). 
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These results are plotted in Fig. 12 as a function of the sensitivity measure in Eq. (69), S di3k = kd/k c . 
Observe how the transition from the classical perturbation to the modified perturbation result, that is, from 
weak to strong localization, occurs at about the sensitivity value 5 mu i ~ 30. Also note that the Monte 
Carlo results agree well with the perturbation results in the limits of strong and weak sensitivity. Figure 13 
illustrates the validity of the perturbation approximations throughout the passbands. In Fig. 13(a) an excellent 
agreement between the simulation results and the classical perturbation solution is observed. Figure 13(b) 
demonstrates the quality of the modified perturbation solution in the passband of a very sensitive system, 
while the classical result grossly overpredicts the localization factor. 

Blade Mistiming: 

The perturbation analyses yield the following localization factors: 

Weak Sensitivity 
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where the plus and the minus signs denote the first and the second passband, respectively. 
Strong Sensitivity 
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with the same sign convention as above. 

Figure 14 illustrates the transition from the classical perturbation approximation of the mid-band 
localization factor to the one obtained by the modified perturbation approach. The sensitivity measure, 
(Eq- (73)), is affected by all parameters, m, k c and k d as well as the passband number. Figure 14 
contains three overlapping Monte Carlo simulation curves, each obtained by varying a different system 
parameter, i.e., m in 1 st passband, k d in 2 nd passband, and finally, k c in 1 st passband. These curves overlap 
nearly perfectly (except for a slight discrepancy at very high sensitivity) and they also agree closely with 
the corresponding simulation curve in Fig: 12 for disk stiffness mistiming. This suggests that 5 = /?'(0) is 
indeed highly suitable as a universal measure of sensitivity, at least for mono-coupled systems. 

Figure 1 5 shows the exponential attenuation as a function of frequency in a system that has low sensitivity 
in the second passband but high first passband sensitivity. In the second passband the classical perturbation 
solution provides an excellent prediction of weak localization. In the first passband the modified perturbation 
approach matches the simulated results, thereby confirming the strong localization approximation. 


Coupling Spring Mistiming: 

Weak Sensitivity, from Eq. (60) 
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Figure 14 Monte Carlo simulations ( ) are used to verify Eqs. (84) ( ) and (86) ( ) . 

Multiple Monte Carlo simulations are obtained by varying m in 1 st passband, k d in 2 nd 
passband and k c in 1 st passband. The standard deviation of blade stiffness mistuning is 5%. 
Note the near overlap of the various Monte Carlo results. 




OK 


(a) s mid 2 e - 1-2 (b) ~ 123 (1 st passband) 

Figure 15 Exponential decay, 7, due to 5% blade stiffness mistuning vs. frequency. In stopbands 

Monte Carlo simulations ( ) agree with the propagation constant for the tuned system 

( ) • In (a) the simulations results overlap the classical perturbation result ( ) , 

Eq. (83), in the 2 nd passband. In (b) the simulation agrees with the modified perturbation 
result ( ) , Eq. (85). 


S c 

7c, mid - y (88) 

Strong Sensitivity is never obtained through mistuning of k c . 

In Fig. 16 we observe the excellent agreement between the classical perturbation result and the Monte 
Carlo simulations in the case of spring mistuning. This was evidenced for all other parameter values, which 
confirms the lack of high sensitivity to spring mistuning. It is interesting to note that the localization factor 
vanishes at the left passband edge. This is reasonable since at that frequency no stretching of the coupling 
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(a) k d = 1 . fc c = 0.1 m = 1 . (b) k d = 1 . fcc = 0.1 m = 1 . (detailed view) 

Figure 16 Localization factor as a function of frequency in a system with coupling spring stiffness 

mistiming of standard deviation 5%. Results from Monte Carlo simulations ( ) agree 

with the classical perturbation solution ( ) in the passband and with the tuned system’s 

propagation constant ( ) in the stopbands. Note the lack of localization at the left 

passband edges and the very weak localization throughout the passbands. 

spring occurs and mistuning of the spring stiffness is inconsequential. 

6.4.3. Summary 

We now sum up the results of this section and attempt to shed light on the meaning of the localization factor. 

The system of Fig. 2 only exhibits low sensitivity to the mistuning of the coupling spring stiffness. This 
spared us the effort of developing a general expression for the localization factor in the case of mistuning of 
parameters connecting two interfaces. Although similar results may be expected for other blade assembly 
models, this should not be generalized to all periodic structures. A counterexample known to the authors 
is a nearly periodic beaded string [18], where mistuning the length of the string segments causes strong 
localization, even though the string connects two interfaces, similarly to the spring in Fig. 2. 

Mistuning of either the disk stiffness or the blade stiffness can lead to strong localization. For disk 
stiffness mistuning we found high sensitivity as k c - 0 or k d - oo. The sensitivity was independent of the 

blade stiffness, k b , the masses m d and m b , and the frequency._ _ 

In the case of blade stiffness mistuning, all parameters fc c , k d and m, as well as the frequency, were 
seen to affect the level of sensitivity. High sensitivity in the first passband was observed as k d oo and as 
m — ► 0 and for either k c - oc or k c - 0. In the second passband the system exhibits high sensitivity as 
m —► oo and k c — > oo, whereas k d by itself cannot directly cause high sensitivity. However, lower values 

of k d cause increased sensitivity in both passbands. 

Figures 12 and 14 depict the transition from the classical to the modified perturbation approximation 
of t in the interval 7 m , d € [0.1, 1]. We have used the subjective terms weak and strong to classify the 
localization in the two limits. According to Figs. 12 and 14 our classification would term 7 < 0.1 as weak 
localization and 7 > 1 as strong localization. In an infinite system with localization the vibration amplitude 
is governed by e~ lN Thus for 7 = 0. 1 the ratio of the amplitudes of adjacent bays is e — 0.90, such 
that on average the energy transmitted from one bay to the next is 82% and 56% is transmitted to the 3 bay. 
For 7 = 1 the average energy transmitted to the next bay is 13.5% and less than 0.25% of the eneigy reaches 
the 3 rd bay! Of course, for the tuned, undamped system 7=0 and 100% of the energy is transmitted. 

The mistuned system whose modes are depicted in Fig. 8 has, at the second midband frequency, a 
sensitivity S bla d e ~ 25. A quick check of Fig. 14 gives 7 m «d ^ 0.2. At 7 = 0.2 the amplitude is decayed 
by 90% by die time it reaches bay number - In |0.1|/0.2 = 11. Of course the vibrations decay in to* 
directions, clockwise and counterclockwise. Modes 14 and 16 in Fig. 8 are closest to the second midband 

and support this. 
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The significance of the value of 7 depends on the total number of bays in the assembly. For example, 
7 = 0.2 indicates 90% attenuation over 1 1 bays and such localization would appear much less significant 
in a 10-bay assembly than it would in a 100-bay assembly. It would thus be useful to define a measure of 
localization that takes the total number of bays into consideration. This can be achieved by introducing the 
localization length or localization scale , N* = I, with which an alternative classification of localization 

could be the ratio In that case an JV-bay assembly would experience 90% localization decay from one 
end to the other when e"* = 0.1, or £ = 2.3, from which the necessary value of 7 (and thus S) could 
be obtained in terms of N. Higher values of jf? would indicate more radical localization. The normalized 
inverse localization length, ^r, could be used as a simple predictive tool. 

7. An Alternative Measure of Sensitivity. 

Results of the proceeding sections indicate that high values of localization are associated with frequency 
ranges in which natural frequencies are closely packed. In Section 4.3 we introduced the concept of modal 
density, the number of natural frequencies per unit frequency. Here we examine the relationship between 
the modal density and localization as a source of an alternative measure of sensitivity. 

Consider the eigenvalue problem: 


{C 0 (u>_,) + 6C(u)j)} Uj = 0, (89) 

where C Q is a circulant matrix describing the dynamics of the full, tuned cyclic assembly and SC is a 
perturbation of that matrix due to mistuning of the system. C 0 and SC are not transfer matrices. It may 
be shown that the first-order perturbation of the eigenvector, Suj, is inversely proportional to the distance 
between the corresponding tuned natural frequency, ljj, and other tuned natural frequencies, u i / j [ 13 ], 

This confirms that if the natural frequencies of a system are closely spaced, the system will be highly sensitive 
to mistuning. 

We attempt to develop a measure of sensitivity based on this information. An expression for the spacing 
of natural frequencies in the cyclic system in Eq. (4) is required. Solving Eq. (30) for n, the natural frequency 
number, and introducing a factor 2 to account for the fact that most natural frequencies are double yields 


n(tD) 
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= 2 


— cos 1 
2tt 



(90) 


where n, a function of ui, is the number of natural frequencies smaller than u>. We define the ratio of the 

derivative of n(Q) to N as the modal density, v (u>), the relative number of modes per unit frequency, and 
obtain: 


d(3 a 

N du 1 -k sin o 


(91) 


The qualitative behavior of Eq. (91) is in some ways similar to that of the localization factor in Eq. (50), e.g., 
near the passband edges (<j — ► 0 or n) the modal density becomes large. However, we were not successful 
in establishing a quantitative relationship between 7 and 77. In the example we examine how correctly 77 
predicts high sensitivity. 
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(a) fc c = 1, 1 st passband (b) m = 1, 1 st passband 

Figure 17 Contours of T) mid , the modal density at midband (Eq. (93)) are plotted in the first passband by 
fixing one parameter and varying the other two. The arrows point in the direction of increasing 


values of 


Example: 

From Eq. (91) we find the modal density of the system in Fig. 2 to be 
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At mid-passband this simplifies to 
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(93) 


Figure 17 displays a parametric study of Eq. (93). Figure 17(b), in which m is held fixed m the first passband, 
confirms qualitative agreement with the equivalent study of the sensitivity measure S mid in Fig. 1 1(c). A 
comparison of Fig. 17(a) and Fig. 1 1(a), in which k c is fixed in the first passband, reveals a vital difference 
between the sensitivity measure, S mid , and the modal density, T, mid . High values of m lead to high modal 
density, which should translate into high sensitivity to mistumng. This high modal density is venfied by 
Fig. A2(a) in Appendix A. Figure 11(a), however, shows low sensitivity values for high values of m The 
Monte Carlo simulations verify that it is indeed the sensitivity measure, 5, that correctly predicts low levels 
of localization in the first passband of systems with high values of m. Although tj correctly yields high 
values of the modal density, it fails to take into account the fact that all the frequencies of the first passband 
go to zero for high values of m, and it is thus a poor predictor of sensitivity to mistumng in this case. We 
conclude that although the modal density is easily calculated and provides useful information regarding 
mistuning effects, it must be used with care. 
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8. Conclusion 

A transfer matrix approach was suggested as an efficient way of modeling blade assemblies in which 
the coupling between remote blades is sufficiently weak to be ignored. We focused on mono-coupled 
assemblies, i.e., assemblies in which only adjacent blade-disk sites are connected and only through one 
degree of freedom. All knowledge about the dynamics of tuned, mono-coupled assemblies is contained in a 
transfer matrix of dimension two, regardless of the actual number of degrees of freedom for each blade-disk 
element. We observed that the frequency domain is divided into passbands, in which waves propagate 
indefinitely, and stopbands, in which attenuated standing waves occur. We reasoned that the number of 
passbands is equal to the number of degrees of freedom for each blade-disk element and that the natural 
frequencies are contained in the passbands. The natural frequencies, mostly double, correspond to a pair of 
complex conjugate, constant interblade phase angle modes. 

In the second part of the paper mistuning was seen to cause a splitting of the double natural frequencies 
and a drastic change in mode shapes. Instead of the extended modes possessed by the tuned system, the 
vibration energy concentrates in relatively few blades with amplitudes far beyond what was predicted by the 
tuned system analysis. These modes of a nearly cyclic system are a manifestation of localization, the effect 
of energy confinement due to scattering of waves at substructure interfaces. 

We suggested a technique by which, based on the model of a single blade-disk element, the sensitivity 
of an assembly to mistuning could be expressed in terms of its parameters, using a sensitivity measure, 5. 
The development of 5 was based on the qualitative changes observed in the system behavior as the expansion 
of the transfer matrix in the small mistuning parameter becomes nonuniform. Perturbation techniques were 
used to generate analytical approximations of the localization factor as a function of 5 in the limits of high 
and low sensitivity. 

To validate the results, Monte Carlo simulations were performed on a popular model of a blade 
assembly. Our conclusion is that our measure of sensitivity, 5, provides an excellent prediction of the levels 
of localization. Not only is the 5 measure simple, intuitive and cost-effective but also amazingly universal, 
i.e., given a statistical distribution in the mistuning parameter, the relationship between 7 m id and S mu j is 
quantitatively the same, regardless of which parameter is mistimed, and is valid for a wide range of the 
system parameters. Approximations of the localization factor showed very good agreement with Monte 
Carlo simulations in the limits of high and low sensitivity, for all frequency values. 
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Appendix A: The Passbands of the 2-DOF Per Bay System: A Parametric Study 

We examine the effect of the parameters rh, kd and k c on the passband/stopband structure of the tuned cyclic 
assembly and hence on the location of its natural frequencies. A physical interpretation of the passband 
edges may be given by recognizing the symmetries exhibited by a cyclic assembly vibrating at the bounding 
frequencies. These symmetries allow us to reduce the problem to that of single blade-disk elements. At the 
edge frequencies adjacent bays are vibrating in phase, a = 0, or out of phase, a = 7r. In the in-phase case, the 
coupling spring does no work and the assembly behaves as if the blades were uncoupled. This model, which 
we shall call the free blade, mimics the dynamics of the left, or lower, passband edges. In the out-of-phase 
case the midpoint of the coupling spring is stationary and the assembly behaves like a series of decoupled 
substructures, each connected to ground through two springs, each of stiffness 2 k c . We refer to this model 
as the fixed blade. It mimics an assembly vibrating at the right, or upper, passband edges. Note that the two 
models hold in both passbands. 



Figure A1 Single blade systems that mimic dynamic behavior at the passband edges. These systems are 
called the free blade system (left) and the fixed blade system (right). 


These two single blade systems are illustrated in Fig. A1 and may be shown to have the natural frequencies 


_2 _ (kd + fh + 1) ± y/ (kd + m + l) 2 — 4m£d 

Wfreei ' 2 " 2rfi 


(Al) 


for the free blade model and 

_ 2 (kd + rh + 4fc c + 1) ± \/ (kd + m + 4 k c + l) 2 — 4 m(fcd + 4 £ c ) 

"*»*.> 2A (A2) 

for the fixed blade model. In both cases the minus sign corresponds to the first passband and the plus sign 
to the second band. 

These single blade systems may be used to examine how the individual parameters affect the pass- 
band/stopband structure of the cyclic assembly. The results are summarized in table Al. The first, and most 
obvious point to be made is that the dimensionless parameter k c has no effect on the lower passband edges. 
In the case of a low k c the upper passband edges approach the lower edges, i.e., the passbands narrow and 
vanish as k c — > 0, as seen in table Al. Less trivial is the case when ic c tends to infinity. The first passband 
motion of the fixed blade system is that of a blade cantilevered at its root with <D = 1 . In the second mode of 
the fixed blade system the disk displacement must be nonzero. Hence as k c — * oo the corresponding natural 
frequency tends to infinity. 

For vanishing disk stiffness k d , the free blade system behaves like a two-mass body free in space. The 
first mode of this system is a rigid body mode and hence the first passband begins at u> = 0. The second 
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natural frequency of the free blade system is u 2 = (1 + fh)/m. The two natural frequencies of the fixed 
system when the disk is not present are 


M3) 


When kd becomes infinite the substructures of the system becomes weakly coupled as evidenced by the 
v anishing of the passbands (see table A 1 ). The spring stiffness k c is obviously immaterial in this case and the 
first natural frequency of both the fixed blade and free blade systems approaches the blade natural frequency, 
tj = l. As in the case of infinite coupling spring stiffness, the high stiffness of the disk degree of freedom 
leads to natural frequencies in the second passband that tend to infinity. 
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Table A1 Passband edges for limiting values of the assembly parameters. The values of (IW, and u>&ee 2 
are given by Eq. (Al) and those of ^ fixed < l d’°) and & Axe d* <r /l > by Eq. (A3). 

Finally, the effect of the mass ratio rh is examined. When m — * 0 the mass of the disk is negligible 
relative that of the blade. The effective inertia of the second mode vanishes with the result that the second 
natur al frequencies of both systems goes to infinity and with it the entire second passband. The first natural 
frequencies of the free and the fixed blade systems approach the distinct values given in Table Al. If, 
conversely, the disk mass is infinitely higher than the blade mass, fh — ♦ oo, the system is equivalent to a 
cantilevered blade and the coupling spring has no effect. The first mode of both systems, in which the masses 
move in phase, becomes a rigid body mode with u> = 0. The second natural frequency of both systems 
approaches the blade natural frequency, 0=1. Thus, both passbands narrow to a single frequency. 

Figure A2 shows a parametric study of the natural frequencies of a 50-blade assembly. Since the 
number of blades is even, the passbands are bounded by the lowest and the highest natural frequencies. 
The most important feature of the system is illustrated in Fig. A2(c). For high values of k c , the natural 
frequencies in the first passband are densely spaced at the right (upper) passband edge — all except the first 
one. This is not surprising since the first mode has all bays vibrating in phase, unaffected by the coupling 
spring, whereas all other modes require some stretching of the very stiff spring. Thus, it is as if the blades 
are mounted on the perimeter of a rigid ring and all frequencies but one approach the natural frequency of 
a cantilevered blade. The first frequency corresponds to a mode in which there is rigid rotation of the ring. 
Note also the implications of this behavior with respect to the midband frequency {a = it /2): The midband 
frequency merges with the right passband edge and is no longer located close to the average frequency of 
the passband. The modal density at the midband frequency tends to infinity. 
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